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Summary 

Numerical investigation has been carried out to evaluate the capability of the Algebraic Reynolds Stress Model 

(ARSM) and the Nonlinear Stress Model (NLSM) to predict strongly curved turbulent flow in a turn-around duct (TAD). 
The ARSM includes the near- wall damping term of pressure-strain correlation (<t> ijtW ), which enables accurate prediction of 

individual Reynolds stress components in wall flows. The TAD mean flow quantities are reasonably well predicted by various 
turbulence models. The ARSM yields better predictions for both the mean flow and the turbulence quantities than the 
NLSM and the k-e (k=turbulent kinetic energy, e=dissipation rate of k) model. The NLSM also shows slight improvement 

over the k-e model. However, all the models fail to capture the recovery of the flow from strong curvature effects. The 
formulation for <J>» w appears to be incorrect near the concave surface. 

The hybrid k-e/ ARSM, Chien’s k-e model and Coakley’s q-co (q = Vk , o) = e/k) model have also been employed 
to compute the aerodynamics and heat transfer of a transonic turbine cascade. The surface pressure distributions and the wake 
profiles are predicted well by all the models. The k-e model and the k-e/ARSM model provide better predictions of heat 
transfer than the q-co model. The k-e/ARSM solutions show significant differences in the predicted skin friction 
coefficients, heat transfer rates and the cascade performance parameters, as compared to the k-e model. The k-e/ARSM model 
appears to capture, qualitatively, the anisotropy associated with by-pass transition. 

Numerical Technique 

An explicit two-dimensional Navier-Stokes solver - RK2D (Kunz & Lakshminarayana, 1992) was used in the 
computation. The RK2D uses a standard 4-stage Runge-Kutta scheme. The fourth-order artificial dissipation is included to 
damp high wave number errors and the second-order artificial dissipation is used to improve the shock capturing. Anisotropic 
scaling of artificial dissipation terms was used. By the use of local velocity scaling, the smoothing was reduced to zero near 
the wall to avoid contamination of the solution by excessive dissipation. Local variable timestepping was also used to 
improve convergence. In the near wall viscous sublayer and buffer layer, the low Reynolds number k-e model is used. The 
matching point for ARSM (or NLSM) and k-e model is at y+= 50 (y+=yux/v, y=distance from the wall, ux=friction 
velocity, v=kinematic viscosity). 
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Resulis-and Discussion 

1) Computation and Analysis of the Flow in a Turn-Around Duct 

Most rocket engine components operate at very high pressures and Reynolds numbers and the flows are fully 
turbulent These flows are often subjected to very large strain-rates such as those arising from the strong streamline curvature. 
There have been many studies on the effects of streamline curvature, however, most deal with mild curvature. Bradshaw (1973) 
reviewed the literature prior to 1973. Recently, Monson et al. (1990) reported detailed measurements in a two-dimensional 
TAD air tunnel at Mach no.=0.1 and Re=10^ or Re=10^. They also calculated this TAD flow using several different versions 
of k-e models and found that only one extended model gave reasonable predictions. The present investigation has been carried 
out with more advanced turbulence models, namely, the ARSM, which was derived from Gibson & Launder (1978) 
Reynolds Stress Model (RSM) by invoking the ARSM assumption, and the NLSM of Shih et al. (1992). A 201x101 
(streamwise x normal direction) H-grid is generated by algebraic method. 

The k-e/ ARSM model was validated against a flat plate boundary layer. The computed values of turbulence 
intensities are in good agreement with Klebanoffs data (1954) for a flat plate boundary layer, as shown in Fig. 1. For the 
TAD flow, all the three turbulence models predict the static pressure coefficients very well on both the inner and the outer 
walls along the bend, as can be seen in the Fig. 2. The pressure loss and static pressure downstream of the bend are also 
predicted reasonably well by all the three models. However, the NLSM predicts a larger separation region on the inner wall 
than the measurement, which can be seen from the underprediction of the pressure coefficient near the exit of the bend. 

Fig. 3(a) shows longitudinal velocity profiles predicted by the k-e model (with Chien’s near wall function slightly 
modified), NLSM and ARSM models. All the models underpredict the velocity near the outer wall. However, the turbulent 
kinetic energy (TKE) and turbulent shear stress (TSS) profiles are predicted reasonably well by all the models. As shown in 
Fig. 3(b) and 3(c), the reduction of turbulent shear stress near the convex (inner) wall is predicted well. The predicted 
enhancement of turbulence near the concave (outer) wall is in agreement with the experiment (Fig. 3(c)). As shown in the 
figures, prediction of turbulence quantities from the ARSM is consistently better than those from the NLSM and k-e models. 

At the exit of the bend (0=180 deg), the data indicate that the flow is separated on the convex wall (Fig. 4(a)). All 
the models predict separation, but underpredict the height of the separation bubble. The prediction of mean velocity by the 
ARSM is slightly better than those from the NLSM and the k-e models. The TKE level near the mid-channel was 

underpredicted significantly by all the models (Fig. 4(b)), this is due to underprediction of the radial component of normal 
stresses near the outer wall. It could be due to the breakdown of the ARSM assumption, i.e., the UjUj/k is constant in the 

flowfield, or the deficiency of the present formulation for the pressure-strain correlation and the dissipation rate of k (i.e., e). 

2) Computation of the transonic turbine cascade flow 

The aerodynamics and blade heat transfer data of the VKI transonic linear turbine guide vane cascade measured by Arts, et 
al. (1990) was selected as the primary test case for present work. The turbine blade shape tested was optimized for a 
downstream isentropic Mach number of 0.9. The downstream isentropic Reynolds number (based on chord length) varied 
from 0.6x 10 6 to 2.1x 10 6 . The total temperature in the free stream was around 410(°k) and the blade surface temperature 
was constant at about 300(°k) for all the test runs. The freestream turbulence intensity varied from 1% to 6%. The inlet flow 
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angle was Pj=0 (deg), chord length C=67.65 (mm), pitch/chord=0.85, stagger angle=55 (deg), axial chord length C ax =37.0 
(mm), design outlet angle p 2 = ^ 4 (deg). For numerical details, see Luo & Lakshminarayana (1993). 

A 129x71 H-grid is used. Blade surface pressure distributions, blade wake, aerodynamic losses and exit flow angles are 
captured very well by all the three models. The blade isentropic Mach no., including the effect of shock wave, is captured 
accurately by all the models, as shown in Fig. 5. Fig. 6 shows the comparison of computed and measured wake profiles. The 
agreement is very good from all the three models. The semi-wake width is shown in Fig. 7. Sudden increase in wake width at 
downstream locations A and B is attributed to shock/wake interaction. It is clear that the correlation shown in Fig. 7 is 
not valid when shock waves are present at the cascade exit. Except for the shock/wake interaction region the wake width 
follows the correlation. In Fig. 8, the predictions of heat transfer by all the models are shown and these are in good agreement 
with the data. The transition is triggered by the shock wave and this has been captured by all the models. The boundary layer 
code (TEXSTAN) underpredicts the heat transfer on the pressure surface and the solution is terminated on the suction surface 
near the inception of separation. Sharma et al. (1982) have observed that most of the disturbance energy is contained in the 
streamwise component of turbulence intensity before transition. During the transition, both the streamwise and the normal 
components grow with the latter component growing at a faster rate than the former, resulting in a decrease of the relative 
magnitude of the streamwise intensity . As can be seen from Fig. 9, this evolution of streamwise and normal turbulence 
intensity within the boundary layer during transition appears to have been simulated qualitatively by the k-e/ARSM model. 
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Fig. 1 Turbulence intensity profiles in the flat-plate 

turbulent boundary layer: experiment by Klebanoff; 
computation by 2-layer k-e/ARSM; 
y=normal distance to the wall, 5= boundary layer 
thickness, Tuj=turbulence intensity components. 
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Fig. 3(a) Longitudinal velocity (U) in turn-around duct, 0=90 
deg.; U m =bulk velocity, 0=angle into bend, y=normal distance 

to the inner wall. 
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Fig. 2 Static pressure coefficient (Cp) on turn-around duct 

inner and outer walls; S=streamwise distance, H=duct 
height, C p =(P-P re f)/(l/2pU^), P re f^static pressure at inlet, 

U m =bulk velocity 
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Fig. 3(b) Turbulent kinetic energy (k) profile 
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Fig. 4(a) Longitudinal velocity in turn-around duct, 0 = 180 deg 
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Fig. 3(c) Turbulent shear stress (<-uv>) profile 
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Fig. 4(b) Turbulent kinetic energy profile 



Fig. 5 Blade isentropic Mach no. distribution for case 

Mur049 (M is> 2=1.02,Rej s> 2=10xl0 6 and Tu„,=l%); 

S = coordinate along blade surface, C = chord length, 
ss = suction surface, ps = pressure surface. 


174 




Fig. 6 Computed and measured wakes at x/C w =1.433 for Fig. 7 Wake width growth downstream of the blade, 
the above case; Por p 02 =total pressure loss, L=semi-wake width, S’/dte=non-dimensional 

x=coordinate along axial chord, y=pitchwise coordinate, distance in exit flow direction 

C^saxial chord length 



Fig. 8 Heat transfer prediction for case Mur222 (M is 2 = 
1.14, Re^ 2=0 55x10 6 and Tu^ =6%); H=heat transfer 


coefficient. 



0.0 0.2 0.4 0.6 0.8 1.0 

Normalized turbulence intensity components 


Fig. 9 The normalized turbulence intensity component 

profiles (Computed by k-e/ARSM) at different axial 

choidwise locations on the suction surface (n is the normal 
distance to the wall and 5 is the boundary layer thickness) 
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